Nature of the Mott transition in Ca 2 Ru04 
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We study the origin of the temperature-induced Mott transition in Ca2Ru04. As a method we 
use the local-density approximation+dynamical mean-field theory. We show the following, (i) The 
Mott transition is driven by the change in structure from long to short c-axis layered perovskite 
(L-Pbca — > S-Pbca); it occurs together with orbital order, which follows, rather than produces, the 
structural transition, (ii) In the metallic L-Pbca phase the orbital polarization is ~ 0. (Hi) In the 
insulating S-Pbca phase the lower energy orbital, ~ xy, is full, (iv) The spin-flip and pair-hopping 
Coulomb terms reduce the effective masses in the metallic phase. Our results indicate that a similar 
scenario applies to Ca2-a;Sra:Ru04 (a; < 0.2). In the metallic x < 0.5 structures electrons are 
progressively transferred to the xz/yz bands with increasing x, however we find no orbital-selective 
Mott transition down to ~ 300 K. 

PACS numbers: 71.27.+a, 71.30.+h, 75.25.Dk, 71.15-m 



The layered perovskite Ca2RuC>4 (Ru 4d 4 , t^e®) un- 
dergoes a paramagnetic metal-paramagnetic insulator 
transition (MIT) at T M it ~ 360 K A similar 

insulator-to-metal transition happens also by application 
of a modest (~ 0.5 GPa) pressure Q and finally when Ca 
is partially substituted by Sr (Ca2- K Sr K Ru04, x < 0.2) 
0, 3- The nature of these transitions, in particular across 
x = 0.2, has been debated for a decade While it 

is clear that a Mott-type mechanism makes the 2/3-filled 
t2 g bands insulating, two opposite scenarios, with differ- 
ent orbital occupations n = (n xyi n xz + n yz ) and polar- 
izations p = n xy — (n xz + n yz )/2, have been suggested. 
In the first, only the xy band becomes metallic, i.e. the 
transition is orbital-selective (OSMT) n and p jump 
from (2, 2) and 1 in the insulator to (1, 3) and —1/2 in the 
metal. In the second, there is a single Mott transition, 




assisted by the crystal-field splitting A — e^^w — e xy > 
13J , similar to the case of 3d 1 perovskites [14j ; p > in 
all phases. To date the issue remains open. Recently, for 
x = 0.2 a novel (xy insulating, n xy = 1.5 and p = 1/4) 
OSMT was inferred from angle-resolved photoemission 
(ARPES) experiments Q, but other ARPES data show 
three metallic bands and no OSMT 0. 

Ca2Ru04 is made of Ru02 layers built up of corner- 
sharing RuOe octhahedra (space group Pbca BE!). 
This structure (Fig. [lj combines a rotation of the oc- 
tahcdra around the c axis with a tilt around the b 
axis. It is similar to that of the tetragonal uncon- 
ventional superconductor Sr 2 Ru04; the corresponding 
pseudo-tetragonal axes x, y and z are shown in Fig. [T] 
The structure of Ca2Ru04 is characterized by a long c 
axis (L-Pbca) above T$ ~ 356 K and by a short one 
(S-Pbca) below Tg. The L- and S-Pbca phases are also 
observed in Ca2- 2; Sr 2 ;Ru04 for all x < 0.2, but Ts de- 
creases with increasing x; for x > 0.2 the system becomes 
tetragonal (for x < 1.5: I4i/acd, c-axis rotations only). 

Because of the layered structure, the ~ xz,yz band- 



FIG. 1: (Color online) Ca 2 Ru0 4 L-Pbca @,[H. The primi- 
tive cell is orthorhombic with 4 formula units; x ~ (a + b)/2, 
y ~ (b — a)/2, z = c are the pseudotetragonal axes. Ru sites 
i at Tj (i — 2, 3,4) are equivalent to site 1 at Ti, with oper- 
ations a — > —a (i = 2), c — > — c (i = 3), b — > b (i = 4) and 
Tj — ► Ti. In the S-Pbca structure the tilting angle is about 
twice as large, while the rotation angle is slightly smaller. 



width, W x 



is about one half of the 
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band- 



width, W xy - Due to the structural distortions, the ti g 
manifold splits into non-degenerate crystal-field states. 
Many-body studies of 3-band Hubbard models show that 
a large difference in bands widths, a crystal-field splitting 
A and a finite Coulomb exchange interaction can affect 
the nature of the Mott transition 0, [l3|> Simple 
models neglect, however, the actual effects of distortions 
on the electronic structure; such effects could be crucial 
14j to the mechanism of the MIT. On the other hand, 
approximate treatments of the many-body effects @ , or 
the neglect of the spin-flip and pair-hopping contribution 
to the Coulomb exchange interaction, could also lead to 
wrong conclusions on the origin of the transition. 
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In this Letter we address the problem by means of the 
LDA+DMFT (local-density approximation + dynami- 
cal mean-field theory) approach (l7j with a continuous- 
time quantum Monte Carlo (QMC) solver 0. This 
method allows us to treat realistically both the material- 
dependence and the many-body effects. We show that 
the Mott transition occurs because of the L — »S-Pbca 
structural phase transition, which is also responsible for 
~ xy orbital order (00). In the metallic phases we find, 
with increasing x, a progressive transfer of electrons from 
xy to xz/yz (p < 0); down to —300 K, we find, however, 
no orbital-selective Mott transition. 

We use the ab-initio downfolding approach based on 
the N-th Order Muffin-Tin Orbital (NMTO) method to 
construct from first-principles material-specific Wannicr 
functions [l4j which span the t2 g bands, and the corre- 
sponding, material-specific, three-band Hubbard models, 
with full local Coulomb interaction [l9| 

H = - X] t mm' C ima C i'm'a + U ^n lm ^n im ± 

+ J2(U -2J - J5„. 
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TABLE I: Ca 2 -xSr !! ,Ru04 (x = 0,0.2): Hopping integrals 
Cm7 me V between sites i' = 1 and i ~ ix + my + nz, 
and (x = 0) crystal-field splitting Ai? ai i/meV=£ Q — ei (a — 
1,2,3) and orbitals at site 1. Orbitals and hopping integrals 
for sites 2, 3 and 4 can be obtained using symmetries (Fig.[T}. 



Here c imcr creates an electron with spin a in the Wannicr 



state m at site i, and n 7 - 



J 



The ab-initio 



parameters t^ m , (Table HJ are the crystal- field splittings 
(i = i') and hopping integrals (i ^ i'), U and J are the 
direct and exchange screened Coulomb interaction. We 
use U = 3.1 eV and J = 0.7 eV, in line with experi- 
mental jUiiil and theoretical [22| estimates. The 
last row in (fTJ) describes the spin-flip and pair-hopping 
Coulomb terms. We solve (J) by DMFT (H] with a weak- 
coupling continuous-time QMC [l8| solver. We retain 
the full self-energy matrix in orbital space, S m m / [lij . 
Our calculations yield the Green function matrix on the 
imaginary axis; we obtain the spectral matrix on the real 
axis by using a stochastic approach [24j]. For limit cases 
(no spin-flip and pair-hopping terms) we perform com- 
parative calculations with an alternative LDA+DMFT 
scheme, based on the projection of LDA Bloch states (ob- 
tained via the projector-augmented wave method [25j . 
Vienna Ab-initio Simulation Package 261 ] ) t o local or- 
bitals and a Hirsch-Fye QMC solver \2^. The pa- 
rameters obtained with the two methods are very similar. 

The nearest- neighbor hopping integrals ££y^cv> ^xltxl^ 
txy^xy progressively decrea se g oing from the ideal tetrag- 
onal Sr 2 Ru04 (see Ref. [29|) to the L-Pbca and then 
the S-Pbca structure of Ca2Ru04; correspondingly, the 
band-width decreases from 2.8 eV (L-Pbca) to 2.5 eV (S- 
Pbca). The cause is the increase in tilting of the RuOe 
octahedra and deformation of the Ca cage, via Ru-0 but 



also Ca-Ru and Ca-0 covalency [14J; for similar reasons 
[bH ]. the crystal-field splitting increases from ^100 meV 
(L-Pbca) to -300 meV (S-Pbca). The crystal-field or- 
bitals are displayed in Fig. [5] For the L-Pbca phase our 
LDA calculations yield p lda ~ 0. This may appear sur- 
prising, since at 2/3 filling due to the difference in band- 
width (W xz / yz ~ 1.5 eV and W xy — 2.8 eV) one might 
expect p < 0. Such effect, however, is cancelled by the 
crystal-field splitting of about 100 meV: since the lowest 
energy state is |1) = 0.932|xy) + 0.259\yz) + 0.255|xz) 
(see Fig. [2]and Table|T|, i.e. close to \xy), neglecting the 
difference in band-width, the crystal-field splitting favors 
p > 0. For the S-Pbca structure we find W xz / yz ~ 1.3 eV 
and W xy — 2.5 eV, i.e. the band-widths decrease by 
about 0.2 eV. The crystal-field splitting is 300 meV, and 
the lowest energy crystal-field orbital, |1), is basically 
identical to \xy). The effect of the crystal-field is stronger 
than in the L-Pbca structure leading to plda ~ 0.37. 

The LDA+DMFT solution of Hamiltonian Q yields 
the following results: For the L-Pbca structure we find a 
metallic solution down to very low temperatures; the or- 
bital polarization is p — at 390 K, i.e. close to the LDA 
value. The self-energy (Fi g. |3l ) exhibits a narrow Fermi- 
liquid regime with kinks [30j; a (lower) estimate of the 
effective masses of the quasi-particles, obtained from the 
slope of Im£ mim (zu; n ) at the first Matsubara frequency, 
is m*/m ~ 5.0 (xy) and m*/m — 4.2 (xz,yz). We find 
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FIG. 2: (Color online) The Mott transition: Correlated bands 
for L- (top) and S-Pbca (bottom), LDA bands (L-Pbca, solid 
lines), and crystal- field orbitals for site 1 (see Table HJ. Pos- 
itive (negative) lobes are red (blue). S-Pbca: |1) is full, \2) 
and 1 3) half-filled; the valence band has character |1). 



similar behavior in Sr2Ru04; in Ca2Ru04 the mass en- 
hancement is larger, because of the narrower band- width. 
Lowering the temperature down to 290 K turns the sys- 
tem into a ferromagnet with (almost) half-metallic be- 
havior and p ~ —0.05 (the occupation of \xy) slightly de- 
creases). The correlated bands for Ca2Ru04 are shown 
in Fig. [5] Along TX we find dispersive bands, while along 
XS the bands become almost flat. 

For the S-Pbca structure the situation is completely 
different. We find a MIT between 390 K and 290 K, 
in very good agreement with the experimental 360 K 
At 580 K the spectral function exhibits a pseudo- 



gap; n xy ~ 1.9 and n x 



s/yz 



1.1, and the polarization is 



already as large as p ~ 0.8. At 290 K the gap is open and 
about 0.2 eV wide (see FigJ3]), while p ~ 1. The most 
occupied state is basically identical to the LDA lowest 
energy crystal- field orbital (Table (JTJ) ) ; the orbital order 
is close to xy ferro-order, with a small antiferro com- 
ponent. LDA+J7 calculations for the anti-ferromagnetic 
phase yield an OO consistent with our results for the 
paramagnetic phase 0, @ . Such orbital order in the S- 
Pbca phase and none in the L-Pbca phase is in line with 
the evolution, across the MIT, of the O Is x-rays absorp- 
tion (XAS) spectra [H, 12| with increasing light incidence 
angle, 9. Our s pec tral functions are also consistent with 
photoemission [l2l [3l| : the quasi- particle peak in the 
high-temperature phase, the Hubbard band at —0.5 eV 



and the multiplcts at -3 eV all correspond to features in 
the experimental spectra. The experimental peak cen- 
tered at ~ — 1.5 eV in [H, 31|, absent in our t 2g spectral 
functions, likely includes a sizable contribution of the O 
bands, which start around those energies. The small gap 
(S-Pbca) is in excellent agreement with electrical resis- 
tivity 0, [H and optical conductivity [l(| data. 

In Fig. [3] we investigate the role played by the spin-flip 
and pair hopping Coulomb terms by comparing results 
with and without those terms. We find that the effec- 
tive masses decrease when the full Coulomb interaction 
is considered; this happens because the average U de- 
creases and the degeneracy of the low energy d 3 and d 4 
multiplets increases. In the insulating phase, such differ- 
ence in multiplet degeneracy appears very clearly as the 
enhancement of the d 4 — > d 3 peak around ~ —0.5 eV. 
The differences in orbital polarization are small. 

Our results show that the Mott transition in Ca2Ru04 
is driven by the L— >-S-Pbca structural transition at T$ 
and occurs simultaneously with orbital order. For the 
L-Pbca structure we find no OO for temperatures well 
below Ts; OO follows, rather than produces, the change 



in structure [32(. These findings are in very good agree- 
ment with resonant x-rays scattering interference data 
which show that the xy OO disappears close to the 
the MIT. They are also in line with the fact (H that a 
tiny pressure (~ 0.5 GPa) is sufficient to make the system 
metallic at 290 K via the S- to L-Pbca transition [2j. 

Can this scenario be extended to Ca 2 - 2; Sr 2 ;Ru04 (x < 
0.2)? Let us examine the limit system x = 0.2, for which 
the L-Pbca phase persists down to 10 K. Neglecting the 



chemical effects of the Ca — > Sr substitution [33j and dis- 
order effects, for the 10 K structure we find crystal-field 
splittings of 81 and 110 meV, slightly (~ 0.1 eV) broader 
t2 g bands than for x = 0, and plda ~ —0.03. The main 
difference with the L-Pbca structure of Ca2Ru04 is in 
the crystal-field orbitals (e.g. (xy\l) ~ (xy\2) <~ 0.66), 
and can be ascribed to the differences in octahedra tilt- 
ing and rotation, and corresponding distortions of the 
cation cage. All these effects stabilize the metallic so- 
lution. With LDA+DMFT (390-290 K) wc find three 
metallic bands, in line with ARPES results from Ref. 0], 
with m*/m ~ 3.7 (xz), 4.4 (yz), 5.6 (xy)] p ~ -0.14. 



While the details slightly differ, depending on x [34(, we 
conclude that for x < 0.2 the temperature-induced Mott 
transition is mostly driven by the change in structure 
L-Pbca -> S-Pbca. We find no OSMT down to 290 K. 

What happens for x > 0.2? For the x = 0.5 structure 
[IH the crystal- field states arc |1) = \xy), and \xz), \yz), 
the crystal-field splitting is small, plda = —0.02, W xy ~ 
2.7 eV, W xz/yz ~ 1.6 eV. LDA+DMFT at 390 K yields 
again a metallic solution with m*/m ~ 4.0 (xz/yz) and 
5.6 (xy), and three metallic bands, in agreement with 
ARPES p 0.15 at 390-290 K. 

Thus, outside the S-Pbca phase we always find a metal, 
in line with transport and optical conductivity data 0, 
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FIG. 3: (Color online) Spectral matrix. Solid lines: Ai,i. 
Dashed lines: ^2,2 and ^3,3. Left: density-density terms only. 
Right: rotationally invariant Coulomb vertex. First row: L- 
Pbca, T=390 K. Second row: S-Pbca, T=290 K. Inset: Real 
(thick lines) and imaginary (thin lines) self-energies. 



EL Ull; with increasing x, n xz + r . 



in line with 

or slightly negative, approaching the 
1/2 of the OSMT scenario Q; we find, however, 



37|; p 



XAS 

P = - 

no OSMT down to 290 K. 

In conclusion, we have studied the origin of the metal- 
insulator transition in Ca2Ru04. We find that it is 
driven by the structural L — ^S-Pbca transition. Two 
mechanisms compete; while a small W xz / yz /W xy band 
width ratio enhances the occupation of the xz jyz orbitals 
(p < 0) and could lead to a orbital-selective Mott transi- 
tion withp = —1/2 0, alarge crystal-field splitting A 3 ,i, 
with |1) ~ \xy) as the lowest e nerg y state, favors xy or- 
bital order and p > 0, as in Ref. [13[ . In the x = L-Pbca 
structure the two effects compensate: Asi/W X y ~ 0.04, 
(l\xy) ~ 0.93 and W xz/yz /W xy ~ 0.54. We find a metal- 
lic solution above and well below Tmit ~ 360 K, with or- 
bital polarization p ~ (no orbital order) at T ~ 390 K. 
At low temperatures the system becomes a ferromag- 
netic metal, in line with (moderate) pressure studies [2j. 
In the x = S-Pbca structure A 3: i/W xy ~ 0.13, size- 
ably larger than for the L-Pbca structure, (l\xy) ~ 0.97, 
and Wxz/yz/Wxy ~ 0.52; the system becomes insulating 
around Tmit and p ~ 1 (~ xy ferro orbital order), in ex- 
cellent agreement with experiments; orbital order follows, 
rather than drives, the transition. Our results indicate 
that this scenario can be extended to Ca2- a; Sr 2 ,Ru04 for 
all x < 0.2. Finally, for the metallic x < 0.5 phases 
we find that, differently than in the crystal-field scenario 
[3], v ~ or negative, slowly approaching the —1/2 of 
Ref. [5j with increasing x, but, down to ~ 300 K, we find 
no orbital-selective Mott transition. 

Calculations were performed on the Jiilich Blue- 
Gene/P. We thank A. Hendricks, T.Koethe, and H. Tjeng 
for sharing with us their data before publication. 
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